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For the laser interferometer space antenna (LISA) to reach it's design sensitivity, the coupling 
of the free running laser frequency noise to the signal readout must be reduced by more than 14 
orders of magnitude. One technique employed to reduce the laser frequency noise will be arm 
locking, where the laser frequency is locked to the LISA arm length. In this paper we detail an 
implementation of arm locking. We investigate orbital effects (changing arm lengths and Doppler 
frequencies), the impact of errors in the Doppler knowledge that can cause pulling of the laser 
frequency, and the noise limit of arm locking. Laser frequency pulling is examined in two regimes: 
at lock acquisition and in steady state. The noise performance of arm locking is calculated with 
the inclusion of the dominant expected noise sources: ultra stable oscillator (clock) noise, spacecraft 
motion, and shot noise. We find that clock noise and spacecraft motion limit the performance of 
dual arm locking in the LISA science band. Studying these issues reveals that although dual arm 
locking [A. Sutton & D. A Shaddock, Phys. Rev. D 78, 082001 (2008).] has advantages over single 
(or common) arm locking in terms of allowing high gain, it has disadvantages in both laser frequency 
pulling and noise performance. We address this by proposing a modification to the dual arm locking 
sensor, a hybrid of common and dual arm locking sensors. This modified dual arm locking sensor has 
the laser frequency pulling characteristics and low-frequency noise coupling of common arm locking, 
but retains the control system advantages of dual arm locking. We present a detailed design of an 
arm locking controller and perform an analysis of the expected performance when used with and 
without laser pre-stabilization. We observe that the sensor phase changes beneficially near unity- 
gain frequencies of the arm-locking controller, allowing a factor of 10 more gain than previously 
believed, without degrading stability. With a time delay error of 3 ns (equivalent of 1 m inter- 
spacecraft ranging error), Time-Delay Interferometry (TDI) is capable of suppressing 300 Hz/x/Hz 
of laser frequency noise to the required level. We show that if no inter-spacecraft laser links fail, arm 
locking alone surpasses this noise performance for the entire mission. If one inter-spacecraft laser 
link fails, arm locking alone will achieve this performance for all but approximately 1 hour per year, 
when the arm length mismatch of the two remaining arms passes through zero. Therefore, the LISA 
sensitivity can be realized with arm locking and TDI only, without any form of pre-stabilization. 



I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) 
is a joint NASA-ESA mission to observe gravitational 
wave signals of astronomical origin. LISA will consist of 
three spacecraft arranged in a triangular formation with 
5 X 10^ m sides (arms). Laser interferometry will be 
used to sense the spacecraft separation to a precision of 
5x = 40 pm/-\/Hz increasing as below 3mHz, giving 
a strain sensitivity of /i ~ 10~^°/VHz. Unlike ground 
based gravitational wave detectors, which are based on 
Michelson interferometers with equal arm lengths, the 
LISA arm lengths vary over time and will be unequal for 
the majority of the mission. Over the period of a year, 
the arm lengths will vary by approximately 3% of the 
total arm length, or by ±76, 500 km [2 . Given this level 
of arm length mismatch, there are stringent laser fre- 
quency noise requirements: coupling of the free running 
frequency noise of a Nd:YAG laser to the signal output 
must be reduced by more than 14 orders of magnitude. 

The frequency noise requirement is expected to be met 
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using three techniques in unison: (1) Pre-stabilization, 
where the laser frequency is locked to either the reso- 
nance of a Fabry-Perot cavity [3J or an mismatched 
arm length Mach-Zehnder interferometer [S]; (2) arm 
locking [3 13 H], a technique based on transferring the 
stability of the LISA arm length to the laser frequency; 
and (3) time delay interferometry (TDI) P [TU] [TT]. a 
post processing technique that synthesizes interferome- 
ters with equal arm lengths by combining phase mea- 
surements with appropriate delays. 

The initial arm locking proposal by Sheard et. al [B] 
showed the round trip propagation delay, r « 33 seconds, 
of the LISA measurement scheme didn't necessarily limit 
the arm locking control system to a low bandwidth, low 
gain system. This is because instantaneous information 
of the laser phase noise is obtained via the local oscilla- 
tor field used in each one-way phase measurement. By 
careful design of the controller, the arm locking control 
system can encompass many interferometer nulls to de- 
liver a high bandwidth, high gain system over the LISA 
measurement band. 

High-bandwidth arm locking has been demonstrated 
in hardware [T^l [131 [H] and analyzed theoretically [S]. 
Dual arm locking [15] , built on the proposal of Enhanced 
Arm Locking by Herz |16j . uses combinations of phase 
measurements from two arms to increase the frequency 
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of the first null of the sensor from 1/r « 30 mHz to 
1/(2At) > 2 Hz , where At is the 1/2 the difference 
in light travel round trip times of the two arms used 
(we define the average round trip time as f). Moving 
the first null to outside the LISA band allows a more 
aggressive controller design below 2 Hz and eliminates 
from the LISA science band noise amplification due to 
the nulls. 

Although there have been many studies of arm locking, 
there has been little analysis of how it would be imple- 
mented in LISA, the operation of it, the performance 
limitations due to noise sources, or the effects associated 
with the changing arm lengths and Doppler shifts. This 
paper is intended to provide a detailed study of these is- 
sues and to predict the performance limitations of arm 
locking in LISA. The analysis is performed for a continu- 
ous system, though the implementation on LISA will be 
digital; we expect the difference in performance will be 
minor. The following issues are analyzed: 

• The arm locking measurement architecture. That 
is, which phasemeter signals should be used in the 
arm locking sensor. 

• Laser frequency pulling due to an error in the 
Doppler frequency estimate. The phase measure- 
ments used in arm locking require an estimate of 
the Doppler frequency to be subtracted. Esti- 
mation error causes the laser frequency to ramp, 
pulled from the starting frequency. Laser frequency 
pulling occurs both at lock acquisition and through- 
out the orbit. To prevent significant laser fre- 
quency pulling, either the Doppler frequency esti- 
mate must be updated regularly or low frequency 
filtering should be implemented. 

• The significant noise sources in arm locking. These 
are: 

— Noise of ultra-stable oscillators (USO's), or 
clocks as we call them henceforth, which en- 
ters at each phasemeter measurement. 

— Spacecraft motion, which is significantly 
larger than the proof mass motion and is 
present in most of the phase measurements. 

— Shot noise, due to the quantized nature of the 
electromagnetic field. 

These effects are investigated using the common and 
dual arm- locking sensors '15], and are written in a gen- 
eral formalism applicable to any sensor. We introduce 
a modification to the dual arm locking sensor that com- 
bines the reduced noise and reduced frequency-pulling of 
common arm locking and the higher gain of dual arm- 
locking. We present a controller design that maximizes 
frequency noise suppression, including a factor of 10 more 
gain than previously assumed |15j , while retaining a con- 
servative phase margin. This controller suppresses laser 
frequency to a level low enough that arm- locking and TDI 



alone, without frequency prestabilization, is adequate for 
the ultimate sensitivity requirement. 

The paper is laid out as follows: In section |IT] the no- 
tation used throughout the paper is introduced, and the 
various possible arm locking sensors previously published 
are listed. The detailed study begins in section [TTT| by ex- 
amining which phasemeter signals should be used in the 
arm locking sensor. In section [TV] laser frequency pulling 
is examined at lock acquisition and in steady state. The 
largest noise sources of arm locking are examined in sec- 
tion [V] and written into a general formalism for the noise 
floor of arm locking in section |VT] In section [VTI] the new 
dual arm locking sensor is presented. In section [VIII| we 
present a design of the arm locking controller. Finally, 
in section |IX| we discuss the performance limitations of 
arm locking to give a realistic estimate of arm locking 
performance in LISA. 



II. NOTATION AND SPECTRA 

In this section we introduce notation, write down the 
arm locking sensors, and the laser frequency noise spectra 
we shall assume for the rest of the paper. 

A. Notation 

Figure [T|shows a simplified schematic of the LISA mea- 
surement scheme and the signal routing for arm locking. 
For simplicity, only two of the three spacecraft are drawn 
and only laser frequency noise is considered. In section VI 
shot noise, clock noise, and spacecraft motion are added 
to the calculation. We label the three spacecraft 1,2, and 
3, and take spacecraft 1 to be the central spacecraft. The 
phasemeter output on each spacecraft is given by propa- 
gating the laser noise source to the phasemeter. To start 
with, assume that all control loops are open. In the fre- 
quency domain, the phase measured at the output of the 
phasemeter on spacecraft 3, facing spacecraft 1 (repre- 
sented by the red circle labelled ^31) is 



(1) 



where 0Lj(w) is the laser phase noise of the laser on 
the jth spacecraft pre-arm locking expressed in units of 
cycles/VHz, and T13 is the light propagation time be- 
tween spacecraft 1 and 3. For simplicity we assume^ 
Tij — Tji. If the laser on spacecraft 3 is phaselocked to 
the incoming light, the closed loop phase noise at the 
output of laser 3 (represented by the red circle labelled 
03) is 



^03(w) 



G3(0.) 



1 + G3(w) 1 + G3(t^) 



(/'Li(^)e 



,(2) 



^ In general, the light propagation time to and from each spacecraft 
will not be equal due to the motion of the constellation. 
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FIG. 1: Schematic of arm locking control loop. Laser fre- 
quency noise is represented by (j}Lj, with j representing the 
number of the local spacecraft. 



where G3{uj) is the frequency response of the phase lock- 
ing controller on spacecraft 3. With the laser on the 
spacecraft 3 phaselocked to the incoming light, but the 
arm locking control loop open {Gi{uj) = 0), the phase 
at the phasemeter output labeled A13 on spacecraft 1 is 
given by 

1 + Gsicj) 

where Pi^^u) is the frequency response of the 1-3 arm, 
including the prompt and delayed signals. The calcula- 
tion of phase 4>Ai2i^) follows a similar argument. The 
frequency responses of the two arms are 



1 



G2{i^) 

1 + G2(w)' 



-22a;Ti3 



(4) 



with G2{i^) the controller gain of the phase locking loop 
on spacecraft 2. In the high gain transponder limit^ 
Gj(w)/(1 -I- Gj{uj)) 1 and 1/(1 + Gj{u;)) for 
j — 2,3, and the frequency responses of the two arms 
are approximately 

Fi3H~l-e-^2c.xi3^ Pi2M«l-e-'2---. (5) 
and the phases that enter the arm locking sensor are 

0A13([j) « (j)Ll{oj)Pl3{uj), (6) 
0A12(w) « 0Ll(w)Pl2(w)- (7) 

We adopt the approximations of (|5]), ([6|, and ([7| in the 
sequel. 



B. Arm locking sensors 

The phase measurements that enter the arm locking 
sensor on the central spacecraft are 



4>A12{^^) 



(8) 



We are interested in the performance of both common 
and dual arm locking. This is because the modified dual 



arm locking sensor, introduced in section VII is a hybrid 
of the two. In common arm locking, phase measurements 
from both arms are simply added. The mapping of the 
two input signals for common arm locking is given by the 
vector 



[1 1] 



(9) 



The open loop noise at the output of the arm locking 
sensor (point Bl in figure [ij is simply 



(10) 



We will be interested in the frequency noise at the laser 
output with the arm locking control loop closed. This 
is the noise at the point Ol in figure [T] With the arm 
locking control loop closed {Gi{uj) ^ 0) the frequency 
noise at the laser output is given by 



<t>oi\ + i^) 



Gl(^)P+(^)0Ll(c^) 



1-f Gi(^)P+(^) 



l + Gi{u;)P+iu;)-' 



(12) 



where P+{uj) = Pi^^uj) + Pi2(w), is the frequency re- 
sponse of the common arm locking sensor, also given in 
table B 3 . 

For dual arm locking, the signal mapping vector is |15j 



1-^ 1 



iuj/S.T 



(13) 



The open loop noise at the output of the arm locking 
sensor (point i31 in figure [TJ is 



and the closed loop output is 

0Ll(t^) 



'/'oib(t^) 



1 + Gi{u)Pd{^Y 



(14) 



(15) 



with Po{lj) the frequency response of the dual arm lock- 
ing sensor. The signal mapping vectors of, and frequency 
responses for, single, common, and dual arm locking are 



^ The phaselocking control loops are expected to have a unity gain 
frequency > 30 kHz, yielding high gain across the LISA science 
band. 



^ In defining the sensors in this paper we have assumed the high 
gain transponder limit, thus wo use the approximation of Pi2(aj) 



and Pi2(w) given in equation [5] 
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contained in table |T] The parameter E{u;) in the dual 
arm locking signal mapping vector is a filter designed to 
ensure stability when the common and difference sensors 
are combined 15j. The values used in this paper are 
given in table |lT] 

The sensors in this paper that contain the addition of 
the two arms (common, dual, and modified dual) have 
an frequency response magnitude up to twice that of the 
single arm locking sensor. To take account of this fac- 
tor or 2 difference, we define the controller for any such 
sensor as Gi(w) = For dual arm locking, which 

has a frequency response of magnitude 2 at / < 1/2At, 
this means that we can design a controller Gl{uj) as if we 
had a sensor with unity magnitude, making the design 
process more intuitive. 



C. Laser frequency noise spectra 

The level and shape of the laser frequency noise spec- 
trum is required for a number of calculations in this pa- 
per. We use three different initial laser noise levels: free 
running laser noise, and laser noise predicted for two 
types of pre-stabilization, Fabry-Perot cavity stabiliza- 
tion and Mach-Zehnder stabilization. The correspond- 
ing laser frequency noise levels are (square root of power 
spectral density): 



Az^fr(/) = 30, 000 X 1 Hz// 

Ai'Fp(/) 



Hz 



(16) 



Hz 



800 X (1 + (2.8 mHz//)2) (17) 

VHz 

Hz 



30 X (1 + (2.8 mRz/fY) 



(18) 



The free running laser noise level quoted here is con- 
servative. A measurement of the beat between two free 
running laser Nd:YAG NPRO lasers is shown in figure |2] 
These results were taken by interfering two free running 
lasers on a beamsplitter (both Lightwave model 126), de- 
tecting one output of the beamsplitter on a photodetec- 
tor and measuring the resultant phase of the beatnote 
using the LISA phasemeter f30]. This measurement is 
always under 25,000 x 1 Hz// Hz/\/Hz which, assum- 
ing the two lasers have identical noise properties, indi- 
cates a free running laser noise of one laser of less than 
18,000 X 1 Hz// Hz/THz. 



III. MEASUREMENT ARCHITECTURE FOR 
ARM LOCKING 
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Frequency [Hz] 

FIG. 2: A measurement of the frequency noise between two 
NPRO Nd:YAG lasers (solid trace), and the assumed free run- 
ning laser noise assumed in this paper (dashed trace). The 
phase of the beatnote was measured using the LISA phaseme- 
ter. The dashed trace is given by equation |16| 



measurements with delays determined by inter-spacecraft 
ranging |17j . This process requires five phase measure- 
ments on each optical bench (with two optical benches 
per spacecraft). These phase measurements are [18]: (1) 
inter-spacecraft measurement, (2) the backlink measure- 
ment, (3) the proof mass to optical bench measurement, 
(4) and (5), the beatnotes of the upper-upper and lower- 
lower clock sidebands. Time delay interferometry forms 
the displacement measurement using the strap-down ar- 
chitecture [19_ , which combines the inter-spacecraft mea- 
surements, the proof mass to optical bench, and the back- 
link measurement to remove spacecraft motion. Clock 
noise in the measurement can be removed using the beat- 
note of clock sidebands [Ml HI HH HI HI • 

Unlike TDI, arm locking requires high bandwidth sig- 
nals (^ 20kHz), in real time, and has significantly less 
stringent noise requirements than the LISA science mea- 
surement. Given the relaxed noise requirements, we 
assume that arm locking will operate with the most 
simple measurement architecture - using only the inter- 
spacecraft phase measurements. A consequence of choos- 
ing this measurement architecture is that both clock noise 
and spacecraft motion will be present in the phase mea- 
surements used for arm locking. In section |V] we see that 
clock noise and spacecraft motion represent the largest 
noise sources for arm locking in the science band. 



The LISA science measurement is performed between 
proof masses on the local spacecraft and distant space- 
craft using the post processing technique, TDI. Time de- 
lay interferometry will be implemented on Earth by form- 
ing linear combinations of low bandwidth (^^3 Hz) phase 



IV. LIMITING LASER FREQUENCY PULLING 

The heterodyne frequency of the inter-spacecraft mea- 
surement is set by the Doppler shift of the laser light. 
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TABLE I: The signal mapping vector and frequency response of different arm locking configurations. Here Tij is the one way 
light travel time of the ijth arm, r is the average round trip time of the two arms, and E{uj) is an filter use d to co mbin e the 
common and difference sensors, given in table |TT] The parameters H+{lj) and H-{lj) are defined in equations 70 and 71 



Configuration 


Signal Mapping 


Frequency Response 


Single 


Ss = 


1 




Ps{lj) = 2i sin (nju) e'''^^'^ 


Common 


s+ = 


1 1 




P+{iu) = 2(1 - cos {Atuj) e"*"") 


Difference 


s_ = 


1 - 




P- (lu) = -2j sin (Atlo) e"^"^ 


Dual 


Sd = 


I _ g(^) I 1 1 


Po{cu) = p+H-§^^P-H 


Modified dual 


Sm ~ 






PmH = P+{lo)H+{lj) - P^{lo)H^{lo) 



TABLE 11: Parameters of the filter contained in the dual arm 
locking sensor. The filter is defined as E{u!) = E\E2- 



Filter 


zeros (radians/s) 


Poles (radians/s) 


Gain 


Si 

E2 


22 = 27r X 10/ Ar 


pi = 27r X 0.1394/Ar 
P2 = 27r X 5/(2Ar) 


91 = Pi 

92 = P2/Z2 



due to the relative motion of the spacecraft. The relative 
velocities in each arm will be up to 18 m/s, correspond- 
ing to a Doppler shift of up to 18 MHz [27]. For arm 
locking to operate stably, this round trip Doppler fre- 
quency must be estimated and subtracted in the phase 
measurements used in the arm locking sensor |26j . In 
the limit of a high gain DC coupled arm locking control 
system, an error in the estimated Doppler frequency is 
compensated for by changing the local laser frequency 
to maintain the desired beat note frequency. In single 
arm locking, this frequency change will appear on the 
light returning from the distant spacecraft 33 s later, ne- 
cessitating a further change by the local laser frequency 
to maintain the desired beat note frequency. The closed 
loop master laser frequency, z/cl, will be changed by the 
error in the Doppler frequency, ^de, each round trip, or 
an average rate of 



St 



2L 



[Hz/s] 



(19) 



For example, if the Doppler frequency can be estimated 
to 100 kHz, the laser frequency will be forced to change by 
1 GHz in 4 days. Such large pulling of the laser frequency 
is undesirable, as it could drive the master laser through 
a mode-hop region, compromising instrument sensitivity. 
The other lasers in the constellation are also at risk of 
being pulled into a mode-hop region as they will be locked 
to the master laser frequency. Additionally, ramping of 
laser frequency combined with scattered light sources can 
couple noise into the science band [2H]- 

We study frequency pulling in two regimes: 1) in 
steady state operation, and 2) at lock acquisition. At 
lock acquisition, the laser frequency can be pulled sig- 
nificantly by an error in the initial Doppler frequency 
estimate and also in the time derivatives of the Doppler 



frequency. 

The solution to this problem explored in detail here is 
to add high pass filters to the arm locking control loop 
at a frequency below the LISA science band. This low 
frequency modification to the controller limits the laser 
frequency pulling at lock acquisition and allows indefinite 
operation of the arm locking control system with an ac- 
ceptable level of pulling. An alternate and closely related 
solution presented by Gath [55], is based on a DC cou- 
pled controller with an additional control loop operating 
at low frequencies to hmit the amplitude of the controller 
signal at these frequencies. This active solution may have 
precision advantages in implementation. 

The arm locking control system operates as follows: 
before arm locking is engaged, measurements of the 
Doppler frequency and the Doppler rate (the first time 
derivative of the Doppler frequency) will be made (Ap- 
pendix E]) and subtracted from the phasemeter measure- 
ment. After the control loop is closed the error in the 
Doppler frequency measurement will cause the laser fre- 
quency to ramp at a rate proportional to the product of 
the error and the step response of the controller. Whilst 
locked, the Doppler frequency estimate will not need 
to be updated. The arm locking control loop will be 
unlocked and re-locked periodically to perform mission 
tasks, such as to change the heterodyne frequencies [?7] . 
At these times the Doppler frequency and its time deriva- 
tives will be known very accurately (as many weeks or 
months of data can be averaged to measure it) and the 
impulse to laser frequency will be much smaller than in 
the first time arm locking is engaged. 

In the following sections we describe the extent of the 
frequency pulling and how to limit it. We shall examine 
the lock acquisition and steady state operation regimes 
separately. We start by calculating the frequency re- 
sponse of Doppler frequency error into the stabilized laser 
frequency. These results are considered in the controller 



design process presented in section VIII where we show 
that with the frequency pulling limited to an acceptable 
level, high gain can be achieved across the LISA science 
band. 
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FIG. 3: Block diagram of the single arm locking control sys- 
tem showing where the Doppler frequency error, vue enters 
the control loop. Ps(tj) and Gi{u}) are frequency responses 
of the sensor and controller, respectively. The laser frequency 
with the control loop open and closed are vol and fch- 



A. Frequency response to Doppler frequency error 

1. Single arm locking 

The block diagram in figure |3]shows where the Doppler 
frequency error enters the single arm locking control loop. 
Here i^ql is the initial (open loop) laser frequency, vcl 
the stabilized (closed loop) laser frequency, and Ps(^) 
and Gi(ijj) are the frequency responses of the single 
arm sensor and controller, respectively. The Doppler 
frequency error enters the phase measurement in the 
phasemeter, where an estimate of the Doppler frequency 
is subtracted from the phase measurement. The closed 
loop frequency will be pulled by 



(20) 



The frequency response of the Doppler frequency error 
to the stabilized laser frequency for single arm locking is 
simply 



-Gi(^) 



l + Giiu;)Ps{Lu) 



(21) 



The frequency pulling due to error in the Doppler fre- 
quency can be examined in the extreme cases of low and 
high loop gain. At high gain. 



Ysito) 



1 



Ps(^) 



(for Gi(t^)P5(w) » 1), (22) 



which, as Ps{o-!) ~ it^r at low frequencies (/ ^ 1 mHz), 
gives the constant-rate frequency pulling described in 
equation [T9j At low gain, 

Ys{lo) ^ G,{lo) (for Gi(c^)P5(^)«1). (23) 

This equation shows that for an AC coupled control loop 
the Doppler error is multiplied by the low-frequency gain. 



2. Common and dual arm locking 

When phase measurements from two arms are used to 
create the arm locking sensor the analysis becomes more 
complex. Figure |4] shows the block diagram for a two 
signal arm locking control loop, such as used in com- 
mon or dual arm locking. Common arm locking would 
use only the signal paths with solid lines, whereas dual 
arm locking uses both solid and dashed lines. There are 
two errors in the Doppler frequency, one enters at each 
phasemeter: fDEi2 and t'oEis- The way the Doppler fre- 
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FIG. 4: Block diagram of the common or dual arm locking 
control system. Common arm locking would use only the sig- 
nal paths with solid lines, whereas dual arm locking uses both 
solid and dashed lines. The Doppler frequency errors that 
enter at the two phasemeters are labelled udeu and i^deis. 
Pi2(tj) and P\3{u}) are the single arm frequency responses. 



quency errors enter the control loop with these sensors 
is less intuitive than in the single arm locking case. The 
Doppler frequency error is added after the individual arm 
frequency responses and before the combination of the 
phase measurements to form the common or dual arm 
locking sensors. 

For common arm locking, the closed loop laser fre- 
quency driven by Doppler frequency error is 



-Gi(cj)(;^DEl2 + i^DEia) 

l + Gi(c^)(Pi2(c^) + Pi3(c^))' 

-Gi(c^)S+Vde 

l + Gi(^)P+(^)' 



where 



Vde = 



^'DE12 J^DEia 



(24) 
(25) 

(26) 



For common arm locking, the frequency pulling is caused 
by the sum of the Doppler frequency errors, t'DE+ = 
'^DEi2 + i^DEi3- The relevant frequency response is 



i^DE+ 



-Gi(a;) 
- Gi{u)P+{Loy 



(27) 
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The closed loop laser frequency driven by Doppler fre- 
quency error in dual arm locking can be given in a similar 
form to equation |25] by replacing the signal mapping vec- 
tor and frequency response 



(28) 



There are two relevant frequency responses to Doppler 
frequency errors in dual arm locking: due the sum and 
the difference in the Doppler frequency errors 



D 



c'fDE+ 
dvcL\D 



I + Gi{uo)Pd{uj) 
-Gi{uj)E{uj) 



(29) 
,(30) 



Si^DE- iujl^T{l + Gi{uj)PD{Lo)y 

where j^de- = '^dei2 ~ '^deis- Due to the additional 
filtering in the difference path of the dual arm locking 
sensor, is inversely proportional to tjAr, causing 

it to dominate the frequency pulling in dual arm locking. 

B. Doppler frequency error 

We now consider what the Doppler frequency error will 
be. There will be an error associated with the initial 
estimate of the Doppler frequency. Also, because the 
Doppler frequency estimate will not be updated whilst 
arm locking is operating, (to prevent adding noise in 
the science band) there will be an error due to the 
Doppler frequency changing throughout the year. Thus, 
the Doppler frequency error must be written as a func- 
tion of time. The Doppler frequency error in the ij arm 
in the time domain is 

i^DEjj(i) = voij + / ^ij{t)dt+ I / aij{t)dt'dt + 
Jo Jo Jo 

(Higher order terms), (31) 

where voij is the initial error in Doppler frequency, jij (t) 
is the error in the Doppler rate (the first time derivative 
of the Doppler frequency), and ctijit) is the error in the 
second time derivative the Doppler frequency. We neglect 
higher order time derivatives of the Doppler frequency 
because they are sufficiently small. 



C. Doppler frequency error at lock acquisition 

The analysis in this section is restricted to the study 
of transients which occur at lock acquisition. We shall 
design the controller such that the transients will occur 
over a period of a few days. Over this period the terms 
7ij {t) , aij {t) will change little and for simplicity we shall 
approximate these terms as constants equal to their ini- 
tial errors. The Doppler error as a function of time is 
then 



t'DEij(t) ~ l^Otj + lOtjt + 



aoijt 



(32) 



with the initial errors in the Doppler rate, change in the 
Doppler rate given by 7oij and aoij. 

In the frequency domain, the Doppler error at lock 
acquisition is 



7o»j 
iuj 



OiOij 

2^' 



(33) 



Because we are interested in common and dual arm lock- 
ing configurations we write 



i^DE- = 1^0- 



7o+ 
iuj 
70- 



0^0+ 

ao- 
2t^2- 



(34) 
(35) 

7012 ± 7013 



where Vf)± = foi2 ± t'oiSi sjid similarly 7o± 
and aQ± = aoi2 ± aois- 

At lock acquisition, the closed loop laser frequency 
driven by Doppler error in common arm locking can be 
found using equations |27] and |34| 



-Gi(^) 
1 + Gi{u)P+{lo) 



-CM.. . ..... „.. (-°^ + l7-S^)^W 

The frequency responses to the error in each term in 



equation 36 written separately are 



V+(^) = 



-Gi(c^) 
l + Gi{uj)P+{uoy 

-Gi{uj) 
iuj{l + Gi{Lo)P+{Lo)y 

Giju) 
2Lo^{l + Gi{Lo)P+{uj)y 



(37) 
(38) 
(39) 



These equations show the error in the initial Doppler 
rate, 70+ and change in Doppler rate, ao+ are integrated 
over time and as such they can cause significant pulling. 

To limit pulling due to vo+j the control loop must be 
AC coupled. The lower unity gain frequency can be de- 
termined from the maximum allowable pulling and the 
error in the Doppler frequency. At the unity gain fre- 
quency. 



|Gi(a;ac)/'+(t^ac)| -1. 



(40) 



At frequencies well below l/r the common arm sensor 
can be well approximated by 



/<0.1mHz 



i2ujT. 



(41) 



If the controller gain is flat at low frequencies, the re- 
sponse is rolled off purely by the arm response, and the 
unity gain frequency and the low frequency gain are by. 



\G\M\ = |Gi,| 



1 



(42) 



where we have substituted G\{uj) — 2Gi{lo). This gain 
gives the factor of amplification of the Doppler knowledge 
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error. Equation 42 can be rearranged to obtain the AC 
coupling frequency for a given low frequency gain 



> 



1 



27rf|Gi, 



(43) 



This result sets the lower bound on the unity gain fre- 
quency based on the allowable gain of the Doppler fre- 
quency error. To give maximum gain, we set the lower 
unity gain frequency to as low as possible, for example 
if we allow \Gi^\ = 1000 then /ac > 4.8 x lO^^ Hz. To 
limit pulling due to 70+ and q;o+, further constraints are 
required. 

At lock acquisition, the frequency responses to Doppler 
errors in dual arm locking can be calculated following a 
similar method. Using equations [29] and [34| the frequency 
responses for the common path of dual arm locking are 



-Gi(^) 



-Giju;) 
iw(l-t-Gi(w)PD(a;))' 

2uj^1 + Gi{uj)Pd{uj))' 



(44) 
(45) 
(46) 



found using equations 30 and [35 



The frequency responses for the difference path can be 

(47) 



-Gi{uj)E{lu) 



iwAr(l + Gi(w)PD(w))' 

Gi{uj)E{lu) 
cu^AT{l + Gi{uj)PD{io))' 
Gi{uj)E{uj) 

i2uj^llT{l + Gi{Lj)PD{uj))' 



(48) 
(49) 



A comparison of the frequency responses to common 
and differential Doppler frequency errors in dual arm 
locking again reveals that it is the differential Doppler fre- 
quency errors which will dominate the frequency pulling 
at start up. This is because the differential path has an 
extra zero at DC and is inversely proportional to At. 



D. Laser frequency pulling at lock acquisition 



below the expected TDI capability and to limit laser fre- 
quency pulling. This design has gain less than unity at 
/ < 4.8 /xHz, implemented with a series of 4 high pass 
filters. Further details can be found in appendix [B] Al- 
though this controller was not designed specifically for 
these sensors, the calculation gives both an estimate of 
pulling that can be expected with each sensor, and a com- 
parison of the two sensors. We calculate the frequency 
pulling that would occur with Doppler frequency esti- 
mates with T = 200 s averaging, as part of arm-locking 
initialization. 



1. Common arm locking 

The accuracy of the Doppler frequency estimation 
needed for common arm locking depends on f , the av- 
erage arm length. Estimation errors for two assumptions 
of the laser frequency noise: —free-running lasers and 
Mach-Zehnder prestabilization 



III 



are given in table 

with details in appendix [X] In the case of free running 
laser noise, the estimates of 70+, and q;o+ have an error 
larger than the maximum value determined by the or- 
bital motion ( appendix [C|. Thus we do not use the laser 
measurement of these quantities, rather we assume them 
to be zero. With prestabilization the error in the mea- 
surement of ao+ is larger than its maximum value and 
thus we assume it to be zero. 

Frequency pulling during lock acquisition is deter- 
mined by the step response of the control system. The 
step response follows from the closed-loop control system 
'^cl(s) according to 



^'stop(^) 



-1 / vci.{s) 



(50) 



where s = iuj and is the inverse Laplace transform 
operator. Step responses from the error in the three 
Doppler derivatives are plotted in figure [5] for both free- 
running and prestabilized lasers. The small stairstep pat- 
tern at the start of acquisition has amplitude i/o+ /2 and 
period r. The maximum pulling from the errors in all the 
deriviatives of 1^0+ is 460 and 90 MHz, for the respective 
cases of free-running and prestabilized laser. 



We now have the tools needed to examine the fre- 
quency pulling in steady state and at lock acquisition. 
The magnitude of the pulling at lock acquisition de- 
pends on the sensor used, the controller design, and the 
Doppler frequency estimate. The estimate is made by 
simply observing the beatnote frequency between the in- 
coming and outgoing fields prior to switching arm lock- 
ing on (Appendix |A]) . In this section we determine the 
pulling that would occur for both common and dual arm 
locking at lock acquisition. We shall use the low fre- 
quency part of the controller in section [VIII| which was 
designed to provide sufficient gain over the LISA science 
band to suppress free-running laser frequency noise to 



2. Dual arm locking 

The error in Doppler estimation for dual arm locking 
is computed the same as in common-arm locking, except 
the relevant delay is Ar = 0.123 s, corresponding to the 
arm length mismatch of the 12-13 arm combination at the 
start of the mission. This smaller delay results in smaller 
error i/q- << vo+, but larger pulling. As with common 
arm-locking, when the measurement errors in 79- and 
ao- are larger than the maximum values taken from the 
orbits, they are estimated as zero. The estimation er- 



rors are tabulated in table III and the step responses are 
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FIG. 6: Step response of the different drivers of Doppler fre- 
quency error for dual arm locking with Doppler frequency 
estimate errors with 200 s averaging, made with free running 
laser noise (top plot) and prestabilized laser noise (lower plot) . 



FIG. 5: Step response of the different drivers of Doppler 
frequency error for common arm locking with Doppler fre- 
quency estimate errors with 200 s averaging made with free 
running laser noise (top plot), and prestabilized laser noise 
(lower plot). 



TABLE III: Initial errors in the common and differential com- 
ponents of Doppler frequency, fo+ ,fo- , Doppler rate, 70+ ,70- , 
and second time derivative of Doppler frequency, ao+ ,ao- , for 
T = 200 s averaging. The columns: free running laser and 
Mach-Zehnder refer to the initial laser noise spectrum, given 
by equations 1 16| and |17| respectively. The numbers with * are 
the values of the Doppler rate and change in Doppler rate at 
the nominal start of the mission; the measurement of these 
parameters with 200 s averaging yields errors larger than the 
mean values. For the differential Doppler frequency estimates 
and time derivatives we have taken At — 0.123 s. 



Parameter 


Free running laser 


Mach-Zehnder 


Units 


fo+ 


6.0 X 10^ 


45 


Hz 




2.3 X 10^ 


0.51 


Hz 


70-1- 


-4.3* 


2.2 


Hz/s 


7o- 


1.2* 


0.02 


Hz/s 


a?o+ 


-0.37* 


-0.37* 


/xHz/s^ 


ao- 


1.0* 


1.0* 


AiHz/s^ 



plotted in figure [6] The maximum pulling with estima- 
tion based on free-running lasers is 13 GHz, unsuitably 
large for NdiYAG lasers that have a typical mode-hop 
free range of approximately 10 GHz. With prestabiliza- 
tion, the pulling is dominated by the ap- term, with peak 
to peak value of 250 MHz. 



E. Laser frequency pulling in steady state 

We now consider frequency pulling in steady state, 
long after lock acquisition. In steady state operation the 
pulling will be driven by the change in Doppler frequency 
which occurs as the constellation orbits the Sun. We de- 
termine the pulling in the time domain using a model of 
the Doppler frequencies over the mission shown in ap- 
pendix |C] 



1. Common arm locking 

The laser frequency pulling in steady state for com- 
mon arm locking is given by the convolution of the (t) , 
the inverse Laplace transform of Y-^(lu), and the common 
Doppler frequency, 

i^CL{t)^y+it)^A+{t), (51) 

where * is the convolution operator and A^{t) is the 
common Doppler shift in the two arms used in common 
or dual arm locking shown in figure |35| The resulting 
pulling of the laser is shown in the upper plot of fig- 
ure [Tj It can be seen that the laser frequency pulling 
with the controller designed in section [VHIj is very mod- 
est, less than 8 MHz peak to peak whilst operating in 
steady state. The pulling is independent of the laser fre- 
quency noise as no Doppler frequency estimates are used, 
and will be an insignificant change compared to the laser 
frequency drift over this period. 
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S. Dual arm locking 

The laser frequency pulling in steady state for dual 
arm locking is given by 



(t), (52) 



where A_(t) is the differential Doppler shift of the two 

arms used by dual arm locking and y^'' (t) and y)^\t) are 

the inverse Laplace transforms of yIj^\uj) and \'^), 

respectively. The j/j^ (t) term in equation 
proportional to the arm length mismatch, t 



52 



is inversely 
nich changes 



significantly throughout the year, dominating the pulling 
in dual arm locking. The expected pulling for dual arm 
locking is shown in the lower plot of figure [7] This shows 
unsustainably large pulling of greater than 10 GHz. The 
pulling which occurs when the arm length mismatch ap- 
proaches zero becomes rapid and it is unlikely that the 
lasers could follow this. Also, the time that the frequency 
pulling is unsustainably large would also be a concern. 
For example, when the difference in the 12-13 arm passes 
though zero, the pulling is above 2 GHz for 35 days. If 
all three arms are available, the central spacecraft in the 
dual arm locking sensor can be switched around to pre- 
vent the arm length mismatch going to zero. In this case 
the pulling would be reduced. 
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FIG. 7: The steady state laser frequency pulling of common 
arm locking (upper plot) and dual arm locking (lower plot). 



driving the lasers into a mode-hop region. At lock ac- 
quisition, with a measurement of the Doppler frequency 
with only 200 s of averaging, the laser frequency pulling 
will be less than 460 MHz, if the lasers are free running, 
or less than 90 MHz if using Mach-Zehnder type presta- 
bilization. This will be further reduced if Fabry-Perot 
prestabilization is selected. In steady state, the laser fre- 
quency pulling is expected to be less than 10 MHz peak 
to peak, insignificant compared to the drift of a laser 
frequency over a similar period. 

In contrast, in the form presented in reference [15^, dual 
arm locking would have larger pulling, posing a signifi- 
cant threat of driving the lasers into a mode-hop region, 
and also of coupling in noise from scattered light sources 
into the science band. At lock acquisition, with a mea- 
surement of the Doppler frequency with only 200 s of av- 
eraging, the laser frequency pulling will be approximately 
13 GHz, if the lasers are free running, or approximately 
250 MHz if using prestabiHzation. While 250 MHz of 
pulling at lock acquisition poses only a small threat of 
pulling the lasers into a mode-hop region, the laser fre- 
quency pulling of dual arm locking in steady state is un- 
sustainable. The laser frequency pulling in steady state 
is dependent on the arm length mismatch which passes 
through zero twice per year for each pair of arms. Thus 
the viability of dual arm locking in the form presented 
in reference !15j would depend critically on the availabil- 
ity of all inter-spacecraft laser links. If all laser links 
were available, the pulling could be reduced by switch- 
ing the central spacecraft when the arm length mismatch 
becomes small. 

The laser frequency pulling in dual arm locking can 
be reduced in three ways, 1) reduce the low frequency 
gain of the controller, which sacrifices the frequency noise 
suppression in the science band, 2) Implement a scheme 
to continuously update the Doppler frequency estimate, 
which would have to be done with a view to not add any 
noise in the science band or, the solution we explore here, 
3) create a modified dual arm locking sensor, a hybrid 
common- dual arm locking sensor designed in section [VH] 
This sensor has the frequency pulling of common arm 
locking, but maintains other control system advantages 
of dual arm locking. 



V. NOISE SOURCES IN ARM LOCKING 



F. Discussion of pulling in common and dual arm 
locking 

Whilst maintaining large laser frequency noise suppres- 
sion over the LISA science band (using the controller de- 
signed in section VIII I , the laser frequency pulling caused 



by Doppler frequency errors in common arm locking are 
relatively modest and do not pose a significant threat of 



Significant noise sources in the arm locking control sys- 
tem are clock noise, spacecraft motion, and shot noise. 
These noise sources are common to all arm locking con- 
figurations (e.g. single, common, and dual). We develop 
a general formalism applicable to all arm locking con- 
figurations. First though, an overview and approximate 
amplitude of each of these noise sources is given. 
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A. Clock noise 

The phase of the photodetector signal is determined 
by comparing it to an onboard clock (i.e. USO) signal. 
Consequently, the noise of the clock signal contributes 
an error in the phase measurement that enters at the 
phasemeter. The amplitude of the clock noise added to 
the phase measurement depends on the fractional fre- 
quency noise of the clock and scales linearly with hetero- 
dyne beatnote frequency, A^- (in units of Hz). The clock 
noise added at the ith spacecraft, when measuring the 
incoming light from the jth spacecraft, is given by 



(t>Cij{f) = AyCi(/) cycles/vHz 



(53) 



where Ci(/) is the clock phase noise normalized to a 1 Hz 
clock frequency, given by 



C^{f) 



2vr/ 



:ycles/Hz\/ Hz 



(54) 



and j/(/) is the fractional frequency fluctuations of the 
clock. The amplitude and frequency dependence of y(/) 
is dependent on the type of clock and operational envi- 
ronment, we use y(/) = 2.4 x 10~^^/v7 I/VHz, cor- 
responding to an assumption of the clock stability in 

LISA Bg. 

Clock noise is correlated for measurements on each 
spacecraft. Thus, clock noise that enters phasemeter sig- 
nals on the central spacecraft can add or subtract coher- 
ently with common or dual arm locking. 



B. Spacecraft motion 

To minimize the disturbances of the space environ- 
ment on the measurement the spacecraft will fly drag free 
around the local proof masses |lj . This means the loca- 
tion of the spacecraft relative to the proof mass is sensed 
and controlled to follow the proof mass. Along the axis 
which the science measurement is made, the spacecraft 
follows the proof mass with an error of approximately 
AX = 2.5 X ^1-1- (//0.3Hz)4 nm/VH^ [29]. The inter- 
spacecraft phase measurements will contain this jitter. 
Even with this spacecraft motion, the arms represent an 
excellent frequency reference, as the fractional length sta- 
bility is still small given the very long arm length. Be- 
cause dual arm locking utilizes the fractional stability of 
the arm length mismatch of two arms, the fractional sta- 
bility is degraded accordingly. Although it is possible 
to subtract this motion in real time using the same al- 
gorithm as used for the LISA science measurement, this 
adds complexity and, as we show in section |VI[ the noise 
performance of arm locking is limited by spacecraft mo- 
tion for only short periods, twice per year. 

If the ith spacecraft moves by an amount of AXij(/) in 
the direction of the jth spacecraft, the associated phase 



shift of the laser field traveling from spacecraft i to space- 
craft j is^ 



4>xtj{f) 



A 



cycles / vHz 



(55) 



where A the laser wavelength. Note that in the propaga- 
tion between the ith and jth spacecraft the laser field will 
pick up phase noise due to the motion of both spacecraft. 



C. Shot noise 

Shot noise is a phase error due to the quantization of 
the electromagnetic field. The phase error of shot noise 
is inversely proportional to the square-root of the opti- 
cal power received from the distant spacecraft. Here, we 
model shot noise as phase error added at the photore- 
ciever. Shot noise is given by jl] 



he 1 

2^AR 



1/2 



cycles/ V Hz 



(56) 



where A is the laser wavelength and Pd is the optical 
power received from the distant spacecraft. Shot noise 
at each ij phase measurement is independent and uncor- 
related. 



VI. QUANTIFYING THE NOISE 
PERFORMANCE OF ARM LOCKING 

A. Calculation of noise in arm locking 

In this section we calculate the noise coupling into the 
arm locking control system building on the calculation in 
section [H] The schematic in figure [8] indicates where all 
three noise sources described in section |V] enter the phase 
measurements. The phasemeter output can be written 
down in a similar fashion to section [ll] by propagating 
each noise source to the phasemeter. Again, to start 
with, we assume that all control loops are open. The 
phase measured at the phasemeter on spacecraft 3, facing 
spacecraft 1 (represented by the red circle labelled Ail) 
is given by 

(I^aM = 0L3(w) - ^Li(c^)e-*""^^ - 0xi3(^)e-^""" - 

0X31 (t^) + 0S31 (^) + 0C31 (w) , (57) 

where (pLj is the laser phase noise of the laser on the jth 
spacecraft pre-arm locking. If the laser on spacecraft 3 is 
phaselocked to the incoming light, the closed loop phase 



^ The spacecraft motion on the central spacecraft will be correlated 
in the two arms to some extent. We expect the correlation to be 
small and neglect it here. 
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FIG. 8: Schematic of dual arm locking control loop showing where the different noise sources enter the system. These are laser 
frequency noise, (pLi, clock noise, (pcij, spacecraft motion, (pxij, and shot noise, (j^sij with i and j representing the numbers of 
the local and distant spacecraft, respectively. 



noise at the output of laser 3 (represented by the red 
circle labelled 03) is 



(58) 



1 + G3(c^) 
G3(C^) 

1 + G3(c^) 

0X31 (w) - 0531 (w) - 0C3l(w) 



where G3{uj) is the frequency response of the phase lock- 
ing controller on spacecraft 3. With the laser on the 
spacecraft 3 phaselocked to the incoming light {G3{uj) ^ 
1), but the arm locking control loop open {Gi{lu) = 0), 
the phase at the phasemeter output on spacecraft 1 is 
given by 

0A13(w) ~ 0Li(w)Pi3(w) - 

(1 + e-'^2a>n3) _ 20;,3i(c^)e-^--- + 

0si3(w) + 053i(w)e"'""^-^ + 

0ci3(w) + 0C3i(w)e-^""", (59) 

where successive terms represent: laser phase noise, mo- 
tion of the spacecraft, shot noise, and clock noise. Simi- 
larly, the phase at point A12 is 

0A12(t^) ~ 0Ll(w)Pi2(w) - 

0xi2(^) (1 + e-^2a>n.) _ 20^21 (c^)e-^--- + 
</.si2(w) + 0S2i(w)e-'""^^-H 
0Ci2(w)+0C2i(w)e-^'""^^ (60) 



The phase measurements that enter the arm locking sen- 
sor on the central spacecraft can be rewritten in the com- 
pact form: 



0A13('^) 



= [Ni-KNs + Nc + Nx], (61) 



where is the laser noise sensed at the central space- 
craft photorecievers, Ng is the shot noise, Nc is the clock 
noise, and Nx is the optical bench displacement noise. 
These are given by 



Nr = 



Nr 



0S12(w) -I- 0S2l(w) 
0S13('^) + (t>S3l{^) 



-Pl2(w)0Ll('^) 
^'l3(^)0Ll('^) 

Ai2(Ci(cc>) + C2(c^)) 

-0X12 (^) (1 + e"»2^-i^) - 20x21 (c^) 

-0X13(t^) (1 + e-'2<.^i^) - 20X3l(^) 



(62) 



where we have dropped the time delays in uncorrelated 
noise terms and assumed the heterodyne frequency on 
spacecraft 2 and 3 is equal to the one way Doppler shift. 

The frequency noise at the laser output with the com- 
mon arm locking control loop closed (point Ol in figure [s]) 
is given by 

0oi| + (t^) = 0Li(w) - Gi(u;)0Bi| + (t^), 

0Ll(t^) 
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— — j^^^^— S+ [Ns + Nc + Nx] , 
1 + Gi(w)P+(u;) 

(63) 

where 4>bi\+ — S+^Ai- The control loop acts differently 
on laser frequency noise from other noise sources. In the 
high gain limit, the laser noise is suppressed by the closed 
loop gain, P^{uj)Gi{uj), whereas the other noise sources 
are maximumly imposed on the laser light. 

The closed loop laser frequency noise for other arm 
locking sensors can be calculated in a similar fashion. 
The vector representation and frequency response of dif- 
ferent arm locking sensors are listed in table |T] The closed 
loop noise at the laser output for the dual arm locking is 
given by 

0Olb(w) = (f)Ll{uj) - Gi{uj)(f)Bl\Dii^), 
_ '/'Ll(^) 

l + GiHFi^H 

^^iM^Sd [Ns + Nc + Nx] , 

(64) 

where (j)Bi\D = S_D<i>Ai. 




Frequency[Hz] 



FIG. 9: Noise of common arm locking at the laser output. 
The total comprises clock noise, spacecraft motion, and shot 
noise. This noise budget is plotted for arm length mismatch 
2Ar — 0.51s, corresponding to the maximum arm length mis- 
match. 



B. The noise floor of arm locking 

Armed with equations |63] and |64] the noise floor of 
common and dual arm locking can be found. In this 
section we neglect laser frequency noise, equivalent to 
assuming the controller has infinite gain. In section [Villi 
an arm locking controller is designed and the total noise 



budget, including laser noise, is plotted in section IX 



In the limit of high controller gain, noise sources other 
than laser frequency noise have little dependence on the 
controller gain, depending primarily on the sensor and 
the transfer function into the phase measurement^. 

The noise floors of common and dual arm locking are 
plotted in figure [9] and 10 respectively, in units of fre- 
quency noise. These figures were generated assuming 
high controller gain, the maximum arm length difference 
that occurs in the mission, 2Ar = 0.51 s, and parameters 



listed in table IV The total noise (dashed black curve) 



is a quadrature sum of shot noise (red curve), spacecraft 
motion (green curve), and the clock noise (blue curve). 
For the parameters used here, clock noise represents the 
limiting noise source of dual arm locking below 20 mHz 
and spacecraft motion represents a noise limit at frequen- 
cies above this. Note that clock noise is linearly depen- 
dent on the heterodyne frequency at each phasemeter 



^ The high controller gain assumption mentioned here is valid over 
the LISA science band as the magnitude of the controller gain is 
always greater than 100. If the controller gain is low, the noise 
coupling will be reduced. 



and this noise budget has the worst combination of het- 
erodyne frequencies that can occur for dual arm locking: 
a maximum difference in Doppler shifts between the two 
arms, A13 — A12 — 29 MHz. In the science band, shot 
noise is always smaller than the both clock and space- 
craft motion, though it dominates above band as clock 
noise and spacecraft motion roll off. 

Also plotted in figures |9] and [TO] is the expected TDI 
capability (light green curve). The TDI capability is re- 
laxed below 3 mHz, where acceleration noise of the proof 
mass becomes significant. The frequency noise curve for 
the expected TDI capability is 



Ai^TDi(/) ~ 300 X (1 + (3 mHz//)2) 



Hz 



Hz 



(65) 



This capability assumes that TDI will be limited by the 
time delay error, currently expected to be 3 ns [32]. For 
a 3 ns time delay error, TDI will be capable of suppress- 
ing frequency noise of 300 Hz\/Hz down to the required 
value of residual frequency noise (« 2 pm/ \/Hz per inter- 
spacecraft hnk) ^35J. 

Comparison of figures |9] and [TO] shows the noise floor 
of common arm locking is lower than that of dual arm 
locking over most of the science band, although it peaks 
higher at frequencies corresponding to nulls in the com- 
mon arm locking sensor. The noise floors of both com- 
mon and dual arm locking are substantially lower than 
the expected TDI capability. The noise floor in figure [TO] 
is the optimal noise performance of dual arm locking, as 
it is plotted with the maximum arm length difference. 
While the noise floor of dual arm locking varies with in- 
verse proportionality to the arm length mismatch, the 
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TABLE IV: LISA parameters and amplitude of noise sources 



Parameter 


Symbol 


Value 


Units 


Average arm length 


L 


5 X lO'-* 


m 


Differential arm length 


AL 


< 76, 500, 000 


m 


Laser wavelength 


X 


1064 


nm 


Doppler shift arm 13 


Ai3 


15 


MHz 


Doppler shift arm 12 


Ai2 


-14 


MHz 


Fluctuations of clock 


yit) 


2.4 X 10"^V\/7 


1/VHz 


Shot noise 


4>sij 


10 


/icycles / %/Hz 


Spacecraft motion 


4>xij 


2.5 X ^1 + (//O.SHz)! 


mcycles/\/Hz 


TDI capability 




300 X (1 + (3 mHz//)2) 


Hz/\/Hz 




FIG. 10: Noise of dual arm locking at the laser output. The 
total comprises clock noise, spacecraft motion, and shot noise. 
This noise budget is plotted for arm length mismatch 2Ar = 
0.51s, corresponding to the maximum arm length mismatch. 



noise floor of common arm locking is largely indepen- 
dent of arm length mismatch. Only the magnitude of 
the peaks at frequencies of / = 1/f and integer mul- 
tiples change with differential arm length change. The 
dual arm locking noise floor is plotted with smaller arm 
length mismatches in figure 11 The lower curve (dotted 
black curve) is plotted with 2At — 0.026s, correspond- 
ing to the minimum arm length that arm locking would 
have to operate with if all inter-spacecraft links are avail- 
able and the central spacecraft was switched to give the 
maximum arm length difference |29j . The noise floor is 
a factor of 100 below the TDI capability at the closest 
point showing that if all three inter-spacecraft links are 
available, the noise floor of dual arm locking is sufficient 
for the entire mission. 

In the case where an inter-spacecraft laser link fails. 



FIG. 11: Total noise of dual arm locking at the laser output 
for different arm length mismatches. The lower curve has 
arm length mismatch of 2Ar = 0.026 s, corresponding to 
the minimum arm length mismatch that occurs during the 
mission if any two of three arms are available to choose from. 
The upper curve has arm length mismatch of 2At = 1.5 x 
10"* s, which is the smallest arm length mismatch that the 
total noise is below the TDI capability (green solid curve). 



only two arms will be available, and therefore the cen- 
tral spacecraft cannot be switched when the arm length 
mismatch becomes small. In this case, the arm length dif- 
ference will pass through zero twice per year and at this 
time the noise floor becomes infinite. The upper curve 
(dashed blue) in figure 
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shows the smallest arm length 
mismatch where the total noise fioor is below the TDI 
capability, which occurs for 2Ar = 150 /.ts, equivalent to 
2AL = 44 km. Here, the clock noise, the dominant effect, 
becomes equal to the TDI capability near 3 mHz. The 
spacecraft motion is also less than a factor of 10 away 
from the TDI capability at this arm length mismatch. 
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With arm length mismatches smaller than this, the TDI 
capability cannot be met with dual arm locking. 

The variation of the dual arm locking noise floor due 
to the changing arm length mismatch can be seen in fig- 
ure 1121 This shows the noise sources at 3 mHz over the 
first two years of the LISA mission assuming only two of 
the three LISA arms are available. In this case, the noise 
floor at 3 mHz is below the TDI capability for the vast 
majority of the time and breaches the TDI capability for 
only short periods, twice per year. It also provides some 
indication of how infrequently and how short a duration 
dual arm locking can not meet the TDI capability in case 
of a critical failure of one arm. The noise performance is 
insufficient to meet the TDI capability for approximately 
1 hour, twice per year. 



VII. MODIFICATION OF THE DUAL ARM 
LOCKING SENSOR 

The benefits of the dual arm locking sensor in terms 
of achieving high gain across the science band as well 
as moving the first null in the sensor to above 2 Hz are 
significant, but the laser frequency pulling that occurs 
in steady state is unsustainable. The frequency pulling 
at lock acquisition is also larger than in common arm 
locking. In terms of noise performance, for the majority 
of the LISA orbit the noise fioor of dual arm locking is 
well below the TDI capability. However, when the arm 
length mismatch is small the noise performance of dual 
arm locking is degraded substantially. Here we introduce 
a modification to the dual arm locking sensor to combine 
the benefits of common and dual arm locking. 

The modified dual arm locking sensor comprises the 
common arm sensor at frequencies below the first null 
of the arm response and the dual arm locking sensor at 
frequencies above the first null, with a smooth transi- 
tion between the two. This sensor provides the stability 
and gain advantages of dual arm locking and recovers 
frequency pulling characteristics and the low frequency 
noise fioor of common arm locking. There are no hard- 
ware changes associated with the control system modifi- 
cation. 



A. Design of the modified dual arm loclcing sensor 

The Bode plot of the modified dual arm locking sen- 

(grey curve) along with the 
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sor is plotted in figure 
dual arm locking sensor proposed by Sutton and Shad- 
dock [T5] (blue curve). A block diagram of the modified 
dual arm locking sensor is shown in figure |13| The sensor 
is located on the 'central' spacecraft and uses the usual 
inter-spacecraft phase measurements from two phaseme- 
ters which measure the interspacecraft signal, </>ai3(i^) 
and 0Ai2(w) (assuming spacecraft 1 is the central space- 
craft). These phase measurements are combined to make 
the common and difference sensors, by taking the sum 



and difference, respectively. These signals are then used 
to construct the modified dual sensor. 

The modified dual arm locking sensor is designed so 
that the common arm sensor dominates below the first 
null of the arm (/ < 1/f) and the dual sensor dominates 
above this frequency. The components of the modified 
sensor are plotted in figure |15[ The frequency response 
of the sensor is given by 

Pm(w) = Fc{uj)P+{cu) + FD{uj)PDicj), (66) 

v- ' "-^ V ■' 

Common Part Dual Part 

where the functions Fc{uj) and Fjj^uj) are filters designed 
to smooth the crossover from the common to dual sensors 
given by 

ga9b{s + Zb) 



Fc{lo) 



Fd{oj) = 



s{s + pb) 



9cgd9eS 



with the parameters given in table [V] 



(67) 
(68) 



TABLE V; Parameters of modified dual arm locking filters 



Filter 


zeros (radians/s) 


Poles (radians/s) 


Gain 


Fc{lo) 




Pa = 


ga = if)-' 




26 = 27r X 5/(13r) 


Pb = 2Tv X 5/(2r) 


Qb = Pb/Zb 


FdH 





Pc = 7/(5f) 


gc = i 







Pd = ll/(20f) 


gd = i 







Pe^2n X l/(90r) 


ge = 1 



Equation 66 can be rewritten as a function of the com- 



mon and difference sensors 

Pa/(w) = P+{u)H+{uj) - P_(w)i/_(w), (69) 

with 



H+{iu) - Fc{uo)+Fd{uj), 

ILUiAT 



(70) 
(71) 



For the sensor in figure [l4| Fcitij) has a pole at DC 
along with a lead compensator with a zero at 5/(13f ) Hz 
and a pole at 5/(2f) Hz. Fo{bj) is a high pass filter 
made from four zeros at DC plus poles at 7/(107rf) Hz, 
ll/(207rf) Hz, and two poles at l/(90f) Hz. The fre- 
quency response of the modified dual arm locking sensor 
is similar to the dual arm locking sensor, with an almost 
fiat response below the first null with a magnitude of 2. 



B. Frequency pulling in modified dual arm locking 



Following the formalism set out in section IV the laser 



frequency pulling in modified dual arm locking is 
-Gi(w)SmVde 



\ + Gi{oj)Pm{^) 



(72) 
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FIG. 12: The noise sources of dual arm locking measured at 3 mHz over the first two years of the LISA mission. This plot 
assumes only two of the three LISA arms are available, preventing the central spacecraft from being switched at small arm 
length mismatch. The heterodyne frequency assumed for the clock noise curve here is pessimistic, as we have assumed the 
worst case that occurs in the mission (A12 — A13 = 29 MHz ) for the duration of this plot. The heterodyne frequency over the 
mission lifetime depends on the Doppler shift, which we have neglected for simplicity. 
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FIG. 13: Block diagram of the modified dual arm locking 
sensor. 




Frequency [Hz] 



FIG. 14: Bode plot of the modified dual arm locking sensor 
(grey curve) and the Sutton and Shaddock [T^ (S&S) dual 
arm locking sensor (blue curve). 



where Sm is the signal mapping vector given in table |T] 
Like dual arm locking, there will be pulling due to both 
the common and differential errors in the Doppler fre- 
quency. The frequency responses written in terms of 
equations 70 and 71 are 



vclIm _ ~Gi{lo)H±{lo) 
- Gi{uj)PM{ujy 



1 



(73) 



where the '+' is used for the common path and the '— ' for 
the difference path. Figure[T6]is a block diagram showing 
where the Doppler error enters the modified dual arm 
locking sensor. 



1. Pulling at lock acquisition 

The frequency responses at lock acquisition can be 
found using equation |73] with the appropriate + or — 
sign along with equations [34] and |35| The frequency re- 
sponses are written compactly as 



Vi?(c.) 



-Gi{uj)H±{lo) 
l + Gi{Lo)PM{ioy 
-Gi{uj)H±{u;) 
ic^(1 + Gi(c^)PmM)' 



(74) 
(75) 
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FIG. 15: Common and dual components of the modified dual 
arm locking sensor. The sum of these gives the modified dual 
arm locking sensor. 
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The improved noise performance of modified dual arm 
locking can be seen by comparing the total noise in fig- 
ure [it] to the upper curve in figure [TT] Figure [17] was 
plotted assuming the controller has infinite gain (so laser 
frequency noise is suppressed to zero). The noise floor 
of modified dual arm locking is similar to that of dual 
arm locking above ~10 niHz but asymptotes to the noise 
floor of common arm locking below this. The noise floor 
of modified dual arm locking at this arm length mismatch 
no longer breaches the TDI capability. 
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FIG. 17: Noise at the laser output for modified dual arm 
locking with arm length mismatch of 2At — 1.5 x 10~*s. 



FIG. 16: Block diagram of the modified arm locking control 
system showing where the Doppler frequency errors, t'DEi2, 
i/DEis enter the control loop. 



VIII. DESIGN OF THE ARM LOCKING 
CONTROLLER 



(76) 



again with the ± set to a '+' for the common path and 
the ' — ' for the difference path. 



C. Noise performance 



The arm locking controller is designed for a sensor that 
has a flat frequency response below 2 Hz with a magni- 
tude of 2, and above 2 Hz has nulls in the response and 
a maximum amplitude of 4. Both the dual and modified 
dual arm locking sensors have this frequency response. 
This factor of 2 is normalized out of the controller by 
adding a 1/2 multiplier at the beginning of the controller. 
Although the controller is designed for the maximum arm 
length mismatch, it will operate stably for smaller arm 
length mismatch. 



For modified dual arm locking, the noise at the laser 
output is given by 



Gi(c^) 
l + PA/(a;)Gi(w) 



A. Goal of the arm locking control loop 

The goal is to design a controller such that the arm 
locking control system has sufficient gain across the LISA 
science band to suppress the laser noise to lower than the 
TDI capability. The required suppression of laser noise, 
(77) S'req(/), IS a ratio of the laser noise before arm locking. 
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vohif), and the TDI capability, 



'S'roq(/) 



l^TDl(/) 



(78) 



If arm locking operates without any pre-stabilization, 
z^OL is given by equation [16] If we assume this level of 
laser frequency noise and a time delay error of 3 ns, the 
required suppression is 



'S'rcq(/) 



30/ [//Hz] 



kHz 



300 X (l + (2.8mHz//)2) 



100/ [//Hz] 
1 + (2.8 mHz//) 



2 for 100/zHz < / < IHz. 

(79) 



If pre-stabilization is implemented the required gain 
will be relaxed. 



B. Constraints on the controller design 

Compared to a standard phase locking control loop, 
the arm locking control loop has two additional design 
constraints. These are: 

1. The controller should have appropriate low fre- 
quency filtering to limit the laser frequency pulling. 

2. The controller must allow for the nulls in the sen- 
sor and the additional phase delay associated with 
them. The nulls in the dual arm locking sensor 
occur at frequencies above 2 Hz. 

These additional constraints limit the achievable gain in 
the LISA science band and necessitate careful design to 
ensure loop stability. 

In the following sections we design a controller in two 
parts, treating the low frequency design (at / < 1 mHz) 
and high frequency design (at / > 1 Hz) separately. Fi- 
nally, we combine the parts to complete the controller 
design. 



C. Low frequency filtering 

In section [IV| the laser frequency pulling that arises at 
lock acquisition and in steady state operation was stud- 
ied. At lock acquisition, it was found that low frequency 
filtering is required to limit frequency pulling due to 
imperfect knowledge of the Doppler frequency, Doppler 
rate, and second time derivative of the Doppler frequency. 
The frequency pulling that occurs for modified dual arm 
locking, calculated in section 
that calculated in section ITV 



VII B I is almost identical to 



for common arm locking. 
The lower bound on the lower unity gain frequency for 
common arm locking (equation 43 1 was designed to be 
/ac > 4.8 X 10~^ Hz. If the Doppler frequency estimate 



has an error of !/o+ = 600 kHz (free running laser for an 
averaging time of 200 s) the pulling due to the Doppler 
frequency error will be under 460 MHz, 4.6% of the mode- 
hop- free range^. 

The frequency pulling at lock acquisition due to er- 
ror in the Doppler rate and second time derivative of 
the Doppler frequency can be limited by further low fre- 
quency filtering. We find that three unity gain high pass 
filters with a corner at 0.8 /zHz placed in series with the 
rest of the controller are sufficient to limit this frequency 
pulling. This low frequency part of the controller limits 
the frequency pulling in steady state to less than 10 MHz 
peak to peak for modified dual arm locking. 



D. Rolling up the gain at low frequencies 

To ensure a stable control system the controller-sensor 
system crosses the lower unity gain frequency, /ac, with 
a slope of /. To optimize the gain in the LISA science 
band we consider various slopes that can be used to roll 
the gain up steeply. To design the roll up we set >30 
degree phase margin at the unity gain frequency, and 
aim for the maximum gain at 100/iHz (S'i.oq(100/iHz) = 
10"^ is the requirement). 

Accounting for the ~30 degrees added by the 0.8 /xHz 
high pass filters, a phase margin of 30 degrees and the 
/ slope for crossing the unity gain frequency, leaves an 
additional 30 degrees of phase at 4.8 /xHz which can be 
allocated to the increased roll up of the controller. We 
find that an additional /^ roll up with 5 zeros at 36.6 /iHz 
and 5 poles at 185 /iHz gives a gain at 100 /iHz of «1900 
and meets the phase requirements. 



E. High frequency design of the arm locking 
controller 



The nulls in the dual arm locking sensor mean careful 
design of the controller is needed to ensure loop stabiHty 
and to limit noise amplification. The nulls in the sensor 
have additional phase associated with them and mean 
that there are many unity gain crossings above the first 
null (located at 2 Hz for maximal arm length mismatch) . 

To design the controller slope above 2 Hz we account 
for all phase delays in the control loop, specify a phase 
margin, then allocate remaining phase from 180 degrees 
to slope of the the controller. The delays in the control 
system are summarized in table |VI[ A 30 degree phase 
margin (6'i„argin) has been chosen at frequencies where 
there are nulls in the sensor (/ > 2Hz for the dual arm 
locking sensor). The delays accounted for are: the de- 
lay associated with the arm locking sensor, ^sensor, the 



® The separation of mode-hops in Nd:YAG non-planar ring oscil- 
lator (NPRO) laser are typically ~ 10 GHz. 



19 



delays in the PZT actuator, Tact, the phasemeter pro- 
cessing, Tpni, and in the interaction with the transponder 
spacecraft and the pre-stabilization loop, denoted Ttians 
and Tps, respectively^. The remaining phase from 180 
degrees to be allocated to the controller is 



-^controller 



= 180- 



margin '^sensor 



^^dclay, (80) 



where ^deiay is the total phase due to delays in the control 
loop given by 

fdelay = 360/(ract + Tpm + Ttrans + Tps). (81) 



F. Phase of the sensor at unity gain points 

In arm locking studies published to date [6j [15] the 
phase delay associated with the nulls in the sensor has 
been set to 6'scnsor — — 7r/2. This phase delay was allo- 
cated because it is the maximum possible phase delay of 
the sensor at the nulls |6]. Here, we show that this al- 
lowance is overly conservative when the controller gain is 
low. To determine the sensor phase at frequencies where 
the total open loop gain crosses unity, 0scnsor|uG(^), we 
solve for the phase of the sensor when the loop gain has 
a magnitude of 1, |Gi(cj)| — 1. The open loop frequency 
response of the control system is 



Gl(c^) = Gi(L^)PM(^)e 



— *fdclay 



(82) 



The common arm locking sensor dominates the frequency 
response of modified dual arm locking above the first null 



Pm(^)| 



« Pd(c^) ~ P+{lo), (83) 
= 2(1 -cos(ATt^)e-'"^). (84) 



We are interested in the nulls in the sensor which occur 
when cos (Arw) = 1. Thus let 



Fm(u;)|,>^«2(1- 
Define 6 — tt — i^t/2, giving 

Pm(^)|/>^=2(1 



(85) 



(86) 



Using the identity 1 



„2ie 



— 2 cos 6*6 we see that 9 is 



the phase of the sensor. At unity gain, |Gi(a;)PA/(t^)| — 

\GI{lu)\ ■ 2cOs6'sensor|uG = 1- That gives cos 6'scnsor I UG = 

l/|2G^(w)|, or for the sensor phase at unity gain 

^sensorluG = " arccos (| 2GJ (w) | - ^ ) For \2Gliiu)\ > 1, 

(87) 

Wc have selected the negative solution for 6'scnsor I UG be- 
cause that gives maximum phase delay, which is of inter- 
est for stability considerations. 



^ The phase delay of the prestabihzation control loop can be re- 
moved, as suggested by Sheard et. al |33l . 




2G,Mj= 10 



FIG. 18: Nyquist plot of the sensor for three different con- 
troller gains, \Gl{uj)\. The phase delay associated with the 
sensor corresponds to the angle from the x-axis to the point 
at which the sensor times gain crosses the dashed line which 
has radius of 1. 



At high controller gain, say |2G|(w)| > 100, the point 
on the sensor which crosses unity gain is very close to 
the null. Accordingly, the phase delay of the sensor at 
the unity gain is ^sensor ~ — 7r/2. For a controller gain or 
|2Gi(cj)| — 1 the peak of the sensor response will crosses 
unity gain. The corresponding phase at this point ap- 
proaches zero (^sensor — > 0). The sensor phase at several 
unity gain frequencies is illustrated in the Nyquist plot 
in figure [18) 

With the phase associated with the sensor at unity 
gain and the other delays in table |VI| we find that ap- 
proximately 60 degrees of phase can be assigned to the 
controller, corresponding to a controller slope of /-0-66 
With this slope, a unity gain frequency of up to 14.9 kHz 
can be achieved. A Bode plot of such a controller is shown 
The magnitude is given by \567 /{if)'~''^^\. 



19 



in figure 

The total phase of the control system (black curve) is a 
sum of the system delays (dotted green curve) , the phase 
of the sensor at unity gain (dashed red curve) , and phase 
associated with the controller (blue solid curve). Note 
the behavior of the system delays and the sensor phase at 
unity gain is complementary at high frequencies. This ef- 
fect allows the control bandwidth to be increased and the 
controller slope to be steeper, giving a factor of 10 higher 
gain at 1 Hz than without this effect. If 6'scnsor — — 7r/2 
allocated, the maximal constant controller slope is /""-^s 
and the control bandwidth is 1.2 kHz, giving a gain at 
1 Hz of 61. 
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TABLE VI: Arm locking system delays 



Type of Delay 


Symbol 


Delay 


Notes 


Phase margin 
Arm locking sensor 


^margin 
^sensor | UG 


30° 

-arccos {\2G*i{lo)\-^) 


Control system margin 

Phase of the sensor at unity gain 


Actuator delay 
Phasemeter processing delay 
Transponding S/C phase delay 
Pre-stabilization phase delay 


'^pm 

'^trans 

'7~ps 


5fis 
2/^s 

(27r30kHz)-^s 
(27r30kHz)-^s 


PZT delay (4^s measured) 
DACs have an inherent l^s delay 
30kHz Unity gain frequency assumed" 
30kHz Unity gain frequency assumed 



"By approximating the phase accrued in the transponding phase 
locking system (and that in the pre-stabilization system) by a time 
delay we underestimate the accrued phase at 10 kHz by 0.3 degrees. 
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FIG. 19: Bode plot of the arm locking controller. The lower 
plot shows the total open loop phase of the control system and 
the three components of which it is comprised. These are the 
controller (-59.4 degrees), the system delays (for parameters 
used in table VI I , and the sensor phase at unity gain, given 
by equation |87[ 
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FIG. 20: Block diagram of the arm locking controller. The 
controller is built from five stages: stage 1 consist of three 
very low frequency high pass filters, stage 2, defines the lower 
unity gain frequency; stage 3, rolls up the gain below the LISA 
science band; stage 4 has two poles in parallel to effectively 
transition between stage 3 and stage 5; stage 5 has 9 poles in 
parallel, with gains individually chosen to generate a slope of 
approximately f'"-^^. 



G. The complete frequency response of the 
controller 

The low and high frequency components of the con- 
troller can be combined to obtain the complete frequency 
response. Analytically, the controller can be designed 
with a precise slope at all frequencies. Here, we present 
a controller using only poles and zeros. A block diagram 
of the controller architecture is shown in figure |20| The 
controller consists of five stages. Stage 1 is the very low 
frequency part of the controller and comprises three zero- 
pole pairs to form unity gain high pass filters at 0.8 /iHz 
in series. Stage 2 sets the lower unity gain frequency 
and has a zero at DC and a pole at 210 /iHz. Stage 3 
is a lead stage, which has five zero-pole stages in series 



to roll up the gain steeply between the lower unity gain 
frequency and the low frequency part of the controller. 
Stage 4 consists of two poles in series and provides the 
transition between the low frequency gain and the shal- 
low slope high frequency part of the controller. Stage 5 
is the shallow sloped part of the controller. It consists of 
nine poles in parallel, with the gain for each pole chosen 
to achieve the required slope of approximately 
The frequency response of the controller is given by 

G*(uj) = ( ^l^V X ( X ( Mi±MV X 

\{s+pii){s + pi2) fr[S+P5kJ' 
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FIG. 22: Plot of the disturbance sensitivity, SdCw) and the 
required suppression, l/5'rGq(i^)- Plotted with 2 At = 0.55s. 



FIG. 21: Open loop frequency response, Gl{u}), of modi- 
fied dual arm locking. Also shown in the magnitude plot 
is the gain required to meet the TDI capability, assuming no 
pre-stabilization. In the phase plot, the red curve indicates 
the total phase of the control loop at unity gain (given by 
equation 90 1 . The arm length mismatch is assumed to be 
2At = 0.51s. 



wit h va lues of the zeros, poles, and gains listed in ta- 
ble |vnp. 

The open loop gain of the control system with the mod- 
ified dual arm locking sensor is given by 

Gl(lj) = Gi(w)PM(w)e~*'"("-^'+"p'°+""™+"p=^ (89) 



The Bode plot of Gl{uj) is plotted in figure 21 Note that 
in this plot we have assumed no prestabilization and thus 
set Tps = 0. From the magnitude plot it can be seen that 
the open loop gain is sufficient to meet the required gain, 
Srcqif)- Note that looking at the open loop phase is 
somewhat deceptive. At high frequencies it appears that 
the open loop phase crosses 180 degrees before 1 kHz, 
making the 14.9 kHz bandwidth control loop unstable. 
In fact, the total phase at the unity gain points is always 
greater than -150 degrees which is indicated by the red 
solid curve, which is given by 



180 

9UG = '^GKuj) X h ^scnsorluG + 

TT 



360/(Tact + 

which is valid for / > l/(2Ar). 



(90) 



H. Closed loop performance of the controller 



The closed loop disturbance suppression function is 



shown in figure 22 This function is given by 



5d(w) 



1 



1 + Gl(l^)' 



(91) 



This shows the suppression of the control loop, as well as 
amplification of the noise at the nulls in the sensor and 
above the control loop unity gain frequency. Also plotted 
is the required suppression to meet the TDI capability. 
Note that the amplification at the nulls is always less than 
a factor of 2 except near the final unity gain frequency, 
near 15 kHz, where the amplitude increases to 5. 



IX. EXPECTED PERFORMANCE OF ARM 
LOCKING 

The performance of arm locking can be predicted us- 
ing the frequency pulling analysis of section |IV[ the noise 
analysis of section |VI[ the controller design given in sec- 
tion |VIII and an assumption of the laser frequency noise. 
Here, we calculate noise budgets for three different initial 
laser noise levels: free running laser noise, and laser noise 
predicted for two types of pre-stabilization, Fabry-Perot 
cavity stabilization and Mach-Zehnder stabilization. 



A. Performance assuming free running laser noise 



® An additional low pass filter in series with the controller is gener- 
ally required to roll off the loop gain at the resonance frequency 
of the laser PZT actuator (near 100 kHz) 



The expected frequency pulling which would occur at 
lock acquisition is shown in the upper plot in figure |23[ 
In this case the Doppler frequency errors are the values 
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TABLE VII: Parameters of the arm locking controller 



Stage 


zeros (radians/s) 


Poles (radians/s) 


Gain (radians/radian) 


1 


zi=0 


pi = 27r X 8 X 10"'^ 




2 


22 = 


P2 = 27r X 210 X 10"^ 


52 = 0.95//ac 


3 


23 = 27r X 36.6 x 10"^ 


P3 = 27r X 185 X 10"'* " 


93 = P3/Z3 


4 




P41 = 27r X 3 X 10"^ 
P42 = 27r X 238 x 10"^ 


54 — P41P42 


5 




P51 = 27r X 3 X 10"^ 
P52 = 27r X 3 X 10"^ 
pg3 = 27r X 3 X 10"^ 

P54 — 2-K X 3 

P56 = 27r X 3 X 10^ 
P56 = 27r X 3 X 10^ 
P57 = 27r X 3 X 10^ 
P58 = 27r X 3 X 10** 
P59 = 27r X 3 X 10^ 


551 = 1.3 X 10-=* 
952 = 3.7 X 10-=* 
053 = 4.2 X 10-=* 

554 = 16 X 10-=* 

555 = 30 X 10-=* 

556 = 69 X 10-=* 

557 = 0.11 

558 = 0.33 

559 = 0.70 



" If arm locking is implemented with prostabilization, p:j = 2tt X 
178 X 10-''. 



of vo+ , 7o+ , and ao+ given in the free running laser col- 
umn in table mil These results are near identical to those 
for common arm locking shown in section [IV| because the 
common arm locking sensor dominates the frequency re- 
sponse at low frequencies. Likewise, the result for pulling 
in steady state of modified dual arm locking, shown in fig- 
ure |24[ is similar to that of common arm locking, with 
less than 8 MHz peak-peak pulling. 



Figure 25 shows the noise budget of modified dual 
arm locking with free running laser noise and an arm 
length mismatch of 2Ar = 0.51 s. For this arm length 
mismatch the laser frequency noise is the limiting noise 
source (the system is gain limited) with the other sys- 
tem noise sources well below the laser frequency noise. 
Figure [25] shows that even without any form of laser pre- 
stabilization arm locking will meet the TDI capability 
across the entire LISA science band. At the most sen- 
sitive frequency of LISA, 3 mHz the frequency noise is 
more than a factor of 4 below the TDI capability. If 
there is not a failure of one inter-spacecraft laser link, 
the dual arm locking central spacecraft can be switched 
when the arm length mismatch becomes small, and arm 
locking alone has sufficient performance to meet the TDI 
capability for the mission. 

The performance of arm locking with an inter- 
spacecraft laser link failure is similar to the no failure 
case. For the majority of the mission the system will be 
gain limited as it is in the no link failure scenario, the dif- 
ference is for a few hours per year the system will become 
noise limited, when the arm length mismatch is close to 
zero. The noise performance will be insufficient to meet 
the TDI capability for approximately 30 minutes, twice 
per year. Note that, this noise limited time could be re- 
duced by either fiying a clock with better stability or by 




FIG. 23: The step responses of different drivers of Doppler 
frequency error for modified dual arm locking. The upper, 
middle, and lower plots assume free running laser noise, Mach- 
Zehnder type prestabilization, and Fabry-Perot cavity presta- 
bilzation, respectively, with the Doppler frequency estimates 
averaged for 200 s. 



implementing the clock noise removal algorithm inside 
the arm- locking sensor^. 



^ However, to remove clock noise using this method the clocks at 
the end spacecraft would need to be phase-locked to the central 
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FIG. 24: The steady state laser frequency pulling of modified 
dual arm locking. The pulling due to the common doppler 
shift is the dominant term. 



B. Performance assuming pre-stabilized laser noise 




10"' 10" 
Frequency[Hz] 



The frequency pulling which would occur at lock acqui- 
sition assuming Mach-Zehnder and Fabry-Perot presta- 
bilization is shown in the middle and lower plots in fig- 
ure 23 respectively. The Doppler frequency errors as- 
sumed for z^o+ 7 7o+ 1 ao+ for the former prestabiliza- 
tion are given in the Mach-Zehnder column in table |III| 
and for the latter are 1^0+ = 1.68 Hz, 70+ = 0.08 Hz/s. 
This shows less than 90 MHz peak to peak pulling over a 
period of 10 days for the Mach-Zehnder prestabilization 
or less than 4 MHz for Fabry-Perot prestabilization. 




10 10 
Frequency[Hz] 

FIG. 25: The noise budget of modified dual arm locking with 
arm length mismatch of Ar = 0.51s. The tperformance was 
calculated with free running laser noise as an initial condition. 

In figures [26] and [27| the total noise after arm locking is 
plotted for no pre-stabilization, Fabry-Perot cavity pre- 
stabilization, and Mach-Zehnder pre-stabilization. Fig- 



spacecraft clock in an analogous way to the lasers. 



FIG. 26: Noise floor of modified dual arm locking with dif- 
ferent initial laser frequency noise conditions: free running 
laser noise (blue curve), Mach-Zehnder pre-stabilization (red 
curve), and Fabry-Perot cavity pre-stabilization (black curve). 
Arm length mismatch of Ar = 0.51s. 



ure [26] is plotted for the maximum arm length mis- 
match, while figure [27] is plotted for the minimum arm 
length mismatch, which will be used if there is no failed 
inter-spacecraft laser links. It is clear that either pre- 
stabilization type in combination with arm locking will 
deliver performance several orders of magnitude better 
than the TDI capability. With either pre-stabilization 
system the performance is limited by clock noise and 
spacecraft motion in the region of a few mHz. This is 
the case even with the largest arm length mismatch. 

Again, the performance of arm locking with a single 
link failure is similar to the no failure case with pre- 
stabilization, except at some frequencies the system will 
be noise limited for the whole mission. The ultimate per- 
formance at 3 mHz, for example, will be entirely dictated 
by the noise sources. The noise floor will exceed the TDI 
capability as the arm length mismatch passes through 
zero. The performance is insufficient to meet the TDI 
capability for about 30 minutes, twice per year. 



X. CONCLUSIONS 

We have performed a detailed analysis of the ex- 
pected performance of arm locking in LISA. This anal- 
ysis includes many of the orbital effects expected on 
LISA: changing arm lengths and Doppler frequencies 
as well as laser frequency pulling due to errors in the 
Doppler frequency estimate. It was found that laser fre- 
quency pulling could be limited to an acceptable level 
by high pass filtering the control loop. We also noted 
that not only the Doppler frequency error will cause 
laser frequency pulling at lock acquisition, but because 
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FIG. 27: Noise floor of modifled dual arm locking with dif- 
ferent initial laser frequency noise conditions: free running 
laser noise (blue curve), Mach-Zehnder pre-stabilization (red 
curve), and Fabry-Perot cavity pre-stabilization (black curve). 
Arm length mismatch of Ar = 0.026s. 



the Doppler frequencies continuosly change, the error in 
Doppler rate can too be significant. In terms of the 
magnitude of the puUing, if arm locking is used without 
prestabilization, the Doppler frequency can be estimated 
to 600 kHz over 200 seconds, giving 460 MHz of pulling 
at lock acquisition for the controller designed here. This 
magnitude of pulling would occur only the first time the 
control loop is engaged, as after LISA has been operat- 
ing over longer periods the Doppler frequency estimate 
would be improved. Less pulling at lock acquisition will 
be achieved if either a longer estimate of the Doppler 
frequency was made, or if arm locking was used in com- 
bination with prestabilization. 

The noise analysis presented here included the ex- 
pected dominant noise sources in arm locking; clock 
noise, spacecraft motion, and shot noise. It was found 
that clock noise and spacecraft motion are the dominant 
noise sources in the LISA science band and that they 
will limit the noise performance when the arm length 
mismatch is small. 

We introduced a new sensor design for the dual arm 
locking sensor, which uses a combination of the com- 
mon and dual arm sensor at frequencies below 1/f and 
the dual arm locking sensor frequencies above 1/f. This 
modified dual sensor has the control system advantages 
of dual arm locking with the frequency pulling charach- 



toristics and low frequency noise performance of common 
arm locking. 

We designed an arm locking controller which max- 
imizes gain in the science band, minimizes frequency 
pulling and has a phase margin of greater than 30 de- 
grees to ensure stability. We noted an effect where the 
arm locking sensor phase changes beneficially near the 
unity-gain frequencies of the control system. With this 
effect taken into account, it was found that the control 
bandwidth could be up to 14.9 kHz and still maintain 
a 30 degree phase margin. This control bandwidth is 
10 times higher than if this frequency dependent phase 
effect is not taken into account. 

Using the noise analysis of section IVT] and the controller 



design given in section VIII we calculated noise budgets 
for different initial laser frequency noise levels. The laser 
frequency noise levels we chose correspond to the noise of 
lasers free running, pre-stabilization to a Mach-Zehnder 
interferometer and to a Fabry-Perot cavity. A key result 
is that arm locking alone can meet the TDI capability. 
With pre-stabilization, the TDI capability would be met 
by more than 2 orders of magnitude across the science 
band. 

We analyzed the expected performance of arm locking 
with and without a failure of an inter-spacecraft laser 
link. If there is no failure, arm locking will have sufficient 
noise performance to meet the TDI capability for the 
entire mission, independent of whether pre-stabilization 
is used or not (though the performance is significantly 
better with pre-stabilization). If there is a loss of one 
inter-spacecraft laser link, the noise fioor of arm lock- 
ing will exceed the TDI- capability for approximately 30 
minutes, twice per year, when the arm length mismatch 
is less than 12 km. This would occur with or without 
pre-st abilization . 
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and delayed (and Doppler frequency shifted) laser frequencies. 
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APPENDIX A: DOPPLER FREQUENCY 
ESTIMATION 

The measurement concept for initial estimation of the 
Doppler rate of a given arm is illustrated in Figure |28| 
The transponded, interfered and detected LISA science 
signal (j)M (M for measured) contains the doppler sig- 
nal vm ~ v/\, where v is the relative velocity between 
spacecraft. The variance of the measured frequency, ct^, 
is defined as the integral of the PSD I'm (/): 

r oo 

= / df ul,{f). (Al) 
Jo 

Averaging over a window of duration T introduces the 
filter H{f) = sinc(/r) = sin(7r/r)/(7r/r). The phase 
measurement is 

^M{t) = ^{t)-^{t-T), (A2) 

where 4)(t) is the inherent phase of the laser output and 
r is the round-trip travel time. Differentiating, the cor- 
responding frequency values are 

VM{t)^v{t)-v{t-T). (A3) 

Equation | A 1 1 then becomes 

/"OO 

= / df v\f)mf)\'\L{f)\\ (A4) 
Jo 

with 

L{f) = 1 - exp(27ri/r). (A5) 
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Doppler estimation error 




Frequency noise RPSD 




FIG. 29: Doppler estimation error cr{T) from averaging the 
spectra {>{/) of Figure [30| over period T. 
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a 




KF<yKln)^fFj2>^\n-\\^ + (n + 1)^ - 2 (n3 + 1) 




Ku^log [ll - n-2|ll - n2|i/n^ ([n + l]/\n - 


Kw 


AV/(\/2rn)^l + n - |1 - 7i| 



TABLE VIII: Error in doppler frequency estimation, a, for 
various spectra of frequency noise Round-trip travel 

time = T, averaging time = T, n — T/t. Note the distinction 
between absolute-value brackets | . . . | and grouping brackets 
or (...). 



Using 



\L{f)\^ = 2[1 - cos(27r/T)] = 4sin^(7r/r), (A6) 



d/ i> sinc^(/T)sin (tt/t). (A7) 



We consider — , k — [— 1,— 1,0], corresponding 
respectively to random-walk frequency noise, flicker fre- 
quency noise, and white frequency noise. These shapes 
are typical for the respective sources: free-running lasers, 
USO phase noise, and unequal-arm Michelson interfer- 
ometry limited by white phasemeter noise. Analytic ex- 
pressions for a are in Table [VTlll and representative nu- 
merical values of and the resulting values of a are 
plotted in Figures [30] and [29] respectively. 
The estimation based on USO noise is a new concept that 
requires some explanation. The USO signals are imposed 
as « 10 GHz modulation sidebands on the carrier, for 
clock correction. Just as the carrier is transponded by 
phase-locking the distant laser to the incoming laser, so 
too the USO sideband can be transponded by feeding 
back the sideband/sideband beat to the phase- locking 
input of the USO. Then at the master (local) spacecraft, 
the measured sideband/sideband beat is representative 
of the Doppler shift. The measurement is made insensi- 
tive to laser frequency noise by subtracting the measured 
carrier/carrier phase, leaving noise only from the clock. 
In this manner, the USO can substitute for an optical fre- 
quency reference for the purpose of measuring Doppler 
frequency. 



FIG. 30: Various assumptions of frequency noise as input 
to estimation of the initial Doppler rate v. The uppermost 
curve F represents the noise in a free-running laser, ^(/) = 
kf/fj — 1 10"' Hz^/^; the flat curve represents white 
laser frequency noise, = = 1 x 10* Hz^''^, and the U 
curve is a typical USO spectrum, vu{f) = ki/o/^/J; where 
fc = 2.4 X 10"^^ and i^o is the laser frequency, 2.8 x 10^"* Hz. 



The error in the average of the fe'th time derivative of 
the doppler shift is a generalization of Equation |A1| 



2 



df {2nff''i>l,{f). 



(A8) 



Generalizing Equation |A7| to allow for j stages of aver- 
aging of the fc'th time derivative. 



2 



df (2nf)^'' I'^sinc^-' 



(/T)sin^(^/r) 



(A9) 



For 1/ = f^, Equation |A9| gives finite ajk only for j > 
1 + p + k. Therefore, to estimate the first and second 
derivitaves of the Doppler shift (fc = 1,2, respectively) 
with white frequency noise (p — 0) requires at least two 
and three stages of averaging (j = 2, 3) respectively. Note 
that j stages of averaging requires data of duration jT. 
For the spectrum of free-running lasers (p = —1), the 
minimum number of averages to compute the first and 
second derivatives are j = 1,2, respectively. We have 
computed ajk for k — (0,1,2) and p = (0,1) with up 
to j = 3 stages of averaging, using the Symbolic Math 
package of Matlab 36 . Sample numeric results are shown 
in Figure [31] 



APPENDIX B: APPROXIMATION OF THE 
SENSORS 

In section |IV| we performed a study of the frequency 
pulling caused by the closed arm locking control loop. 
Due to the high computational time required to simulate 
the frequency pulling using the sensors in there exact 
form, the simulations spanning over a period of days af- 
ter lock acquisition were performed using approximations 
of the sensors. These approximations are very good be- 
low the null frequencies, that is for / < 10 mHz in com- 
mon arm locking, and / < 2 Hz for dual arm locking. 
This frequency band is responsible for the most signifi- 
cant pulling occurs. 
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a for 1th derivative of Doppier shift 
v(f) = K/f, K = 30 l<Hz*rthz, tau = 33.00 s 




FIG. 31: Error in the estimate of the first time-derivative 
of the single-arm doppier shift, assuming the noise of free- 
running lasers, figure [30] For averaging times T < r = 33 s, 
there is no advantage to more than one stage of averaging. 
For T > 33 s, two stages of averaging gives less error than 
one, but there is no advantage to three stages. The estimation 
error from two stages, each of duration T = 100 s, is 8kHz/s. 
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■ Dual Sensor 

■ Approximation of Dual Sensor 
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FIG. 33: Dual arm locking sensor and approximation used in 
step responses. 
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■ Common Sensor 

■ Approximation of Common Sensor 
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FIG. 32: Common arm locking sensor and it's approximation 
used in step responses. 



In the study of pulling in common arm locking we ap- 
proximate the frequency response of the sensor by 



4s 



s-f 2/r' 



(Bl) 



which is plotted in figure [32] along with the common arm 
locking sensor. 

In the study of pulling in dual arm locking and modified 
dual arm locking we approximate the frequency response 
of the sensor by 



PuHIa 



(B2) 



which is plotted in figure [33] along with the dual arm 
locking sensor. 

The controllers used in calculating the frequency pulling 
in section |IV| were based on the controller design pre- 
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FIG. 34: Arm length mismatch of the sets of two arms. Cal- 
culations in this paper were performed using arm 12-13. 



sented in section |VIII| The shape of the common arm 
locking sensor differs significantly from that of the dual 
arm locking sensor (for which the controller was de- 
signed). To calculate the frequency pulling in common 
arm locking the controller response was modified to ac- 
count for this. The controller used in common arm lock- 
ing frequency pulling calculations was simply 



Gl{oj) 



approx 

where Gi(a;)* is given by equation 



(B3) 



APPENDIX C: ORBITS 

The orbit data used in this paper was provided by Peter 
Gath [31]. Plots of the relevant parameters are shown in 
figures [34] [35] and [36] 
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FIG. 35: (i) Common Doppler frequency, (ii) common 
Doppler rate, and (iii) 2nd time derivative of common Doppler 
frequency for the different arm pairs. 
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FIG. 36: (i) Difference Doppler frequency, (ii) difference 
Doppler rate, and (iii) 2nd time derivative of difference 
Doppler frequency for the different arm pairs. 



